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For the first time a proper comparison of the average depth of shower maximum (Y max ) published by 
the Pierre Auger and Telescope Array Observatories is presented. The A max distributions measured 
by the Pierre Auger Observatory were fit using simulated events initiated by four primaries (pro¬ 
ton, helium, nitrogen and iron). The primary abundances which best describe the Auger data were 
simulated through the Telescope Array (TA) Middle Drum (MD) fluorescence and surface detector 
array. The simulated events were analyzed by the TA Collaboration using the same procedure as 
applied to their data. The result is a simulated version of the Auger data as it would be observed 
by TA. This analysis allows a direct comparison of the evolution of (Y max ) with energy of both 
data sets. The (Y max ) measured by TA-MD is consistent with a preliminary simulation of the Auger 
data through the TA detector and the average difference between the two data sets was found to be 
(2.9 + 2.7 (stat.) + 18 (syst.)) g/cm 2 . 
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1. Introduction 

Composition is a central key to understand the origin of ultra-high energy cosmic rays. Large 
efforts in developing new detectors and analysis procedures have been made recently in order to 
improve our knowledge about the abundance of particles arriving at Earth. At the highest energies 
(E > 10 18 eV) the depth of shower maximum (A niax ) is one of the most robust observables available 
to infer the composition. Currently, the Pierre Auger and the Telescope Array observatories measure 
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Fig. 1 . (X max ) as measured by the Pierre Auger (left) and Telescope Array (right) Collaboration [2,3]. The 
colored lines denote predictions of air shower simulation (note that different models are shown in the left and 
right panel, only Sibyll2.1 is the same). The black line on the right panel is a straight-line fit to the TA data. 


X max using fluorescence detectors. Despite the use of the same detection principle, a direct compari¬ 
son of the data published by both collaborations is not straightforward. 

The TA Collaboration publishes (X max ) values obtained from distributions folded with detector 
resolution and efficiency. The interpretation of the data is made possible by the publication of the 
Monte-Carlo prediction for proton and iron nuclei also folded with detector resolution and efficiency 
(Fig. 1, right). In the Auger Collaboration only certain shower geometries are selected for sampling 
almost unbiased X max distributions. The corresponding (X max ) values are presented in the left panel 
of Fig. 1. In the Auger analysis, each selected geometry allows a wide enough range of X max values 
to be observed within the fluorescence detector field of view boundaries. We will refer to this event 
selection as ’fiducial selection’. Besides that, the Auger Collaboration published <x(X max ) with detec¬ 
tor resolution unfolded. This procedure allows the interpretation of the data (i.e. (X max ) and cr(X max )) 
using Monte Carlo predictions without the need to fold the detector properties into the predictions. 
The advantage of the TA analysis is that it does not require removing as many events, since this tech¬ 
nique does not apply a fiducial selection. 

The work reported here is a common effort of the Auger and TA Collaborations with the aim 
to provide the cosmic ray community a direct comparison of the (X max ) measurements taking into 
account the different approaches of each collaboration. Indirect comparisons of TA and Auger re¬ 
sults were published in the first report of these series [1]. The disadvantage of indirect comparisons 
is that they depend on the particular hadronic interaction model that is used. The current analysis 
was performed in the following way. The Auger X max distributions were fitted by a combination of 
four primary nuclei (proton, helium, nitrogen, iron) using events from air shower simulations. The 
abundances which best fit the Auger data were simulated through the TA-MD detector and analyzed 
by the TA Collaboration using the same procedure as applied to their data. This procedure resulted 
in the Auger data folded into the TA-MD detector. The Auger (X max ) folded with TA-MD analysis is 
shown in this paper in comparison to the TA-MD data as it is usually published. 
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Fig. 2. Number of selected events for the Auger (solid red line) and TA (blue dashed line) analyses. The ratio 
of events is given in the lower panel. 


2. Data Samples 

The analysis presented here is based on the data measured by the Pierre Auger Observatory in 
the period from the 1st of December 2004 until 31st of December 2012. All measured events were 
analyzed as explained in reference [2]. The events were selected in order to guarantee good mea¬ 
surement conditions and a high-quality reconstruction. After that, the fiducial selection was applied. 
In total 19,947 events were considered for further analysis. The X max values of these events were 
sampled in 18 energy bins starting at log (E/eV) = 17.8. 

From the Telescope Array we use hybrid data collected with the MD fluorescence telescope 
and surface detector array over the period from the 27th of May 2008 to the 2nd of May 2013. The 
reconstruction and analysis applied to the data is described in [3]. The number of events which passed 
all cuts is 438, for which the mean X max is shown in 12 energy bins above log (£/eV) = 18.2. 

The number of events used for this comparison presented here is shown in Fig. 2 and the X max - 
resolution of the two experiments is presented in Fig. 3. As can be seen, the resolutions after cuts are 
comparable but it is worthwhile noting that the resolution quoted for the MD does not contain effects 
from the detector calibration and atmospheric monitoring. The systematic uncertainties on the X max 
scale are compared in the right panel of Fig. 3 and they are < 10 g/cm 2 and 16 g/cm 2 for the Auger 
and TA analyses respectively. 

3. Analysis 

Due to the different analysis approaches of the TA and Pierre Auger Collaborations it is not 
possible to compare the published values of (X max ) and cr(X max ) directly. Whereas the moments of 
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Fig. 3. X max resolution (left) and systematics of the X max scale (right) for the Auger and TA analyses. 



Fig. 4. Moments of the fitted X mm distributions using QGS JetII-03 (black markers) and X mdx moments mea¬ 
sured by the Pierre Auger Collaboration (red circles with statistical error bars), see text. 


the X max distribution published by the Pierre Auger Collaboration are close to the true moments 
(moments of / true in Eq. (1)), the TA collaboration published the (A max > folded with the effects of the 
detector response and reconstruction (moments of f D b s in Eq. (1)). 

The relation between the true and observed X max distribution is 


foUXZ x ) = 


-f 

Jo 


/true(X max ) £(X max ) R(XZ x ~ X max ) dX, 


( 1 ) 


i.e., the true distribution / true is deformed by the detection efficiency e and smeared by the detector 
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resolution R that relates the true X max to the reconstructed one, X£ ax . 

To be able to perform nevertheless a comparison of the two results, we need to establish what 
(Xmax)obs would look like in the TA detector given the X max distribution measured by Auger. For this 
purpose, we convolute a parametric description of / true that is based on the Auger data with the TA 
detector simulation and apply the same reconstruction and analysis chain used for the TA data to this 
simulated data set (see [5] for a previous description of this method). 

Technically, the parametric description of the X max distribution is realized by providing a set 
of composition fractions as a function of energy that describe the X max distributions measured by 
Auger. These fractions arc obtained as described in [4] by a log-likelihood fit of templates of X max 
distributions for different nuclear primaries as predicted by air shower simulations using a particular 
hadronic interaction model. It is worthwhile noting that the detector acceptance and resolution at a 
given primary energy depend mainly on X max itself and only weakly on the primary particle type or 
hadronic interaction model via the invisible energy. Therefore, for the analysis presented here, it is 
only important that the resulting composition mix describes the data well and not which fractions of 
primaries arc needed or which hadronic interaction model is used to obtain the model of the ’true’ 
A niax distribution. 

Here we used QGSJetII- 03 [6] and a mix of four primary particles (proton, helium, nitrogen and 
iron) to obtain a model of true X max distribution based on the Auger data. QGSJetII- 03 is not included 
in the set of models studied by the Pierre Auger Collaboration to infer the primary composition [4] 
because it gives a worse description of LHC data than the re-tuned version QGSJetII- 04 [7]. How¬ 
ever, with neither version of QGSJetII it is possible to find a composition mix that gives a perfect 
description of the X max distributions measured by Auger. The first two moments of the best fits with 
QGSJetII- 03 and the Auger data are shown in Fig. 4. As can be seen, there is a good agreement 
regarding (A max ), but there arc deviations between the fitted and observed width of the distribution. 

Ideally, this analysis should be performed with a combination of composition and hadronic in¬ 
teraction model that fits the Auger data well, such as Sibyll2.1 [8] or Epos-LHC [9] (see discussion 
in [4]). However, due to the lack of large air shower libraries other than QGSJetII- 03 within the TA 
Collaboration, we performed the analysis with this model for practical reasons. Since the deviations 
between the moments of the data and the ones of the fitted distributions arc on average at the 5 g/cm 2 
level, this approach is expected to give only a small bias in the predicted observed distributions. 

In detail, the analysis proceeds as follows: the composition mix is processed using the Telescope 
Array hybrid reconstruction analysis software. Showers arc generated by CORSIKA and the trigger 
response of the surface detector is simulated. The generated longitudinal shower profile is fitted to a 
Gaisser-Hillas function to determine the shower parameters and a profile based on these parameters 
is generated. The TA fluorescence detector response including atmospheric, electronics, and geomet¬ 
rical acceptance is then simulated. Subsequently the event geometry is fitted via the fluorescence 
profile and the shower-detector plane is measured. A fit to hybrid shower geometry is performed 
which combines the timing and geometric center of charge of the surface detector array, with the 
timing and geometry of the fluorescence detector that observed the event. This step is what makes the 
event a hybrid event. If either the surface or fluorescence detector fail to trigger in an event, it is not 
processed any further, otherwise the shower profile is fitted via a reverse Monte Carlo method where 
the atmosphere, electronics, and geometrical acceptance of the shower arc fully simulated. 

The resulting effect of the folding of protons and the parametric Auger distributions with the TA 
detector response, reconstruction and analysis on the (X max ) of Auger is shown in Fig. 5. As can be 
seen, the observed mean is smaller than the unbiased mean. 
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Fig. 5. Effect of the MD detector acceptance on A max . The (X max ) of an X max distribution describing the 
Auger data before and after the MD acceptance are shown as solid squares and circles respectively. The error 
bars denote the statistical uncertainties of the Auger result in case of the squares and the statistical uncertainties 
due to the limited MC statistics in the case of the circles. 


4. Results and Discussion 

The (X max ) as measured by TA using the MD fluorescence telescope and the Auger result folded 
with the TA acceptance are shown in Fig. 6. Their compatibility is quantified with a bin-by-bin com¬ 
parison excluding the highest-energy data point of each experiment which are at different energies. 
Using only the statistical uncertainties yields a^ 2 /Ndf of 10.7/11 with P(y 2 > 10.7|11) = 0.47. The 
average difference of the data points is (2.9 ± 2.7 (stat.) ± 18 (syst.)) g/cnr 2 with a^f 2 /Ndf of 9.5/10 
(P = 0.48). It can be concluded that the two data sets are in excellent agreement, even without ac¬ 
counting for the respective systematic uncertainties on the A max scale. Flowever, in the present study 
we did not take into account a possible difference in the energy scale of the two experiments. The 
comparison of the energy spectra at the ankle region suggests that the energy scale of TA is about 
13% higher than the one of the Pierre Auger Observatory [10]. But since the elongation rate of the 
folded Auger data is small (~ 19 g/cm 2 /decade), the effect of such an energy shift on the comparison 
is expected to be at the level of a few g/cnr 2 . For a more precise evaluation it is required to take into 
account the energy dependence of the acceptance of TA. Nevertheless, it is to be expected that the 
increased difference between the two data sets once the energy scale shift is taken into account will 
be much smaller than the systematic uncertainties on the A max scale of < 10 g/cm 2 and 16 g/cnr 2 for 
the Auger and TA analyses respectively. 

5. Conclusions 

In this paper we presented a comparison between the data on (A max ) as measured by the Pierre 
Auger and Telescope Array Collaborations. An adequate comparison was achieved by taking into 
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Fig. 6. Comparison of (X max ) as measured with the MD of TA (blue squares) and the (X max ) of the Auger data 
folded with the MD acceptance. The data points were slightly shifted horizontally for better visibility. In the 
case of the Auger points (red circles), the inner error bars denote the statistical uncertainty of the measurement 
and the total error bar also includes contributions from the limited statistics of simulated events used for the 
folding. The colored bands show the systematic uncertainties of the X max scales of each experiment. 


account that the (X max ) published by Auger is corrected for detector effects, whereas the (X max ) pub¬ 
lished by TA includes detector effects. In the future, we intend to improve the parametric description 
of the Auger X max distributions and the evaluation of the effect of the relative energy scale uncertainty. 
Nevertheless, from the preliminary comparison presented here we conclude that the data of the two 
observatories are in good agreement. 
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